The effect of a scanning flat fold mirror on a CMB B-mode experiment 
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We investigate the possibility of using a flat-fold beam steering mirror for a CMB B-mode exper- 
iment. An aluminium flat-fold mirror is found to add ^0.075% polarization, which varies in a scan 
synchronous way. Time-domain simulations of a realistic scanning pattern are performed, and the 
effect on the power-spectrum illustrated and a possible method of correction applied. 
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I. INTRODUCTION 

The prevailing ACDM cosmological model has been 
very successful in explaining many observations of the 
universe. However, inflation - exponential expansion 
within the first ~ 10~ 35 s - has many uncertain details; 
indeed it lacks strong confirmation. Inflationary models 
predict a background of gravity waves [l|, H, [2CJ, [H, l23j 
which produce a curl-like, or "B-mode" pattern of po- 
larization on the cosmic microwave background (CMB) 
on scales greater than 1°, which cannot be produced by 
density perturbations [III HH • The amplitude of this sig- 
nal is characterized by the tensor-to-scalar ratio, r, which 
is constrained by current measurements of the tempera- 
ture power spectrum and temperature-polarization cross- 
correlation to a value r < 0.20 JT2| ■ To put tight con- 
straints on inflationary models [2J we need to detect the 
primordial B-mode signal, or constrain it for r < 0.01. 
This requires measurement of polarization fluctuations 
at the level of 10 nK. 

To measure this signal level is challenging, and good 
control of systematics is required [l(J. A number of 
ground-based (QUIET [p, PolarBear QJ]), balloon- 
based (EBEX [7], Spider 4]) and satellite-based (Planck 
[Til li~7j ]) experiments have been, or are close to being, 
fielded with the raw sensitivity to detect the B-mode sig- 
nal if it is present at a reasonable r. These experiments 
are all designed to have low, and stable, instrumental 
polarization (IP), and often include a modulation strat- 
egy to assist with reducing systematics and mitigating 
against 1/f noise. The majority of these experiments 
hold their optics stationary with respect to the receiver, 
and scan the cryostat and optics together; a notable 
exception being the Tenerife experiment [f|. Scanning 
the cryostat places requirements on the stability with re- 
spect to accelerations and tilts. Although cryogenic re- 
ceivers have been demonstrated that are stable at differ- 
ent tilt angles @ , the sensitivity to accelerations has not 
been demonstrated. A system can be imagined, however, 
which combines a stationary receiver with a fold mirror. 
Besides the simplification of the cryogenic requirements, 
the large reduction in the mass of moving components 
would make the design of the steerable mount much sim- 
pler. 

We investigate such a system in this paper. In Sec- 
tion |TT] we describe the system, and review the production 
of polarization from reflections. In Section lflTl we describe 



a fiducal scan strategy and simulation pipeline. Results 
from this pipeline are then presented and discussed, along 
with other potential problems in Section IIVI 



II. METHOD 

For this investigation, we consider a horn fed receiver 
and compact range antenna (CRA), placed in a ground- 
screen, nominally observing the horizon. The output 
beam of the CRA is then steered to the sky with a single, 
flat, fold mirror. Although the flat mirror would be large; 
~ 4 X 4 m, such a system can be considered. The flat 
mirror is then tipped and tilted such that the beam of 
the telescope moves across the sky to allow observations. 

Off- normal reflections, as caused by this folding, induce 
a polarization. The effect is well known and is the prin- 
ciple behind, for example, the "pile-of-plates polarizer" . 
The Fresnel equations give the reflection coefficients for 
incident angle 9i and transmission angle Qt for an inci- 
dent wave with E-field parallel to the boundary, i?y and 
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For a metal plate, we use the relation for the complex 
index of refraction, 



(n + in) 2 = e' + ie'' 



(3) 



For a general dielectric, the conductivity a is given by 
a = uje"eo where uj — 2irf and / is the radiation fre- 
quency. For a good conductor, we can assume e' <C e", 
and so 
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Plugging this into Snell's law allows calculation of R» 
and R±. For a pure aluminium plate, I /a = p = 
28.2 x 10~ 9 tim, we find = 0.9970 and R ± = 0.9985. 
Thus an unpolarized input reflected with 9i = 45 would 
have an overall polarization of 0.075% This ignores the 
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phase-delay from the reflections, which elliptically polar- 
izes the reflected beam, producing both circular polar- 
ization and rotating the plane of polarization. Using a 
Drude-Lorentz model to calculate n and n for aluminium 
at 90GHz 24] , we calculate that the plane of polariza- 
tion is rotated by an additional ~ 10 -9 degrees when 
the phase delay is included. Ignoring this is then valid 
as this is a much smaller than the effect from differen- 
tial reflection (~ 0.2°). The reflection is a total power to 
polarization conversion (T — > Q, U); in the context of a 
CMB polarization measurement the total power is 2.7 K, 
and so the induced polarized signal is then ~200 /iK, 
higher than even the EE polarization spectrum. How- 
ever, the signal is changing consistently with the scan, 
and so is varying on different scales from the CMB sig- 
nals of interest and so could, in principle, be removed. 
In the next section, we describe the pipeline with which 
we investigate the effect on CMB reconstruction and Ci 
estimation. 

We note in passing that it is possible to construct a 
two mirror system that provides a fold with no net po- 
larization. Consider a beam travelling along the x axis. 
The first mirror, mounted at 45° to the yz plane, then 
reflects the beam so it is travelling along the y axis, and 
the second, mounted at 45° to the xy plane, reflects that 
beam along the z. As the reflection axes have switched, 
both input E x and E y receive an attenuating factor of 
RuR±. Any net polarization is then due to differences 
between the conductivity of the mirrors. If the two mir- 
rors are fixed relative to each other, and allowed to rotate 
about the x axis, the beam can be steered. A second set 
of mirrors, in a similar configuration, can be added to al- 
low full beam steering around the sky. This is illustrated 
in figure [1] 



III. SCAN STRATEGY AND PIPELINE 

We simulate observations with the following four step 
pipeline. Firstly, time ordered data (TOD) is simulated 
with 100Hz sampling for four 3-hour observation periods, 
corresponding to a day of observations, and the sky lo- 
cation (RA, Dec), telescope pointing direction (az, el), 
and observed E% and Ey recorded. The input CMB 
maps were produced by simulating a power spectrum 
with C AMB [Tj, and producing a realization with the 
Healpix[Notel| [6j synfast routine. This input all-sky 
map was also convolved with a Gaussian beam with 
FWHM of 5.5 arcmin. Parameters used in the simulation 
are listed in Table [TJ Note that we include the monopole 
in this simulation, with a temperature of 2.7 K and so 
El ~ El ~= 1.3 K. 

It is generally desirable for ground-based and sub- 
orbital experiments to scan the sky at constant eleva- 
tion to minimize the effect of variable airmass. We fol- 
low this and consider the proposed Clover scan strat- 
egy, t he optimization of w hich is fully described else- 
where North and Johnson| . The strategy is to perform 



TABLE I. Simulation parameters. Detector positions are 
given in (az, el) in the telescope frame. 



Parameter 



Value 



Telescope location 
Patch centre 
Patch radius 
Azimuthal scan speed 
Nominal elevations 
Observing frequency 
Data rate 



Atacama desert 
9h00m,-44°00' 
10° 
0.5 deg/s 
45, 60 deg 
90 GHz 
100 Hz 



Resitivity of aluminium 


28.2 x 10" 


9 dm 


Input r 


0.1 




Detector 1 position 


(0,0) arcmin 


Detector 2 position 


(14.8, 14.8) 


arcmin 


Detector 3 position 


(29.6, 29.6) 


arcmin 


Detector 4 position 


(36.9, 7.4) 


arcmin 


Detector 5 position 


(51.7, 22.1) 


arcmin 



azimuthal scans at fixed elevation (nominally 45°) as the 
chosen sky patch rises for around 1 hour, continually ad- 
justing the azimuthal center to achieve as circular a patch 
as possible. Once the patch has risen sufficiently through 
the scan, the elevation is increased, nominally to ~ 60°, 
and the patch continues to rise for another 3 hours. The 
same strategy is used as the patch sets, giving four dis- 
tinct data streams for a single patch during a 8-12-hour 
observation period. 

The second step in the pipeline is to take the pointing 
from the TOD data, determine the angle of the flat fold- 
ing mirror, rotate the observed E% and Ey to align with 
the plane of the mirror, apply the calculated R\\ and R±, 
and then rotate back. The TOD is then naively binned 
to produce T, Q and U maps, following the map-making 
scheme outlined in Q , and then written out as a Healpix 
format FITS files. We make maps at N s ^ c = 512 so that 
the multipole range I ~ 100 to 1000 is retrievable, and 
avoid holes in the maps by simulating a five pixel focal 
plane. The pixel angular positions from boresight are 
listed in Table UJ 

Finally, the angular power spectrum (APS) coefficients 
Ci are estimated using the anafast Healpix routine. We 
apply a circular mask of radius 10° with a simple cosine- 
squared taper of width 2° at the edge. 



A. Determination of angles 

In order to determine the effect for a mirror at an ar- 
bitrary orientation, two angles are required; the angle of 
incidence, and the angle, /3, between the mirror and 
the horizon as seen by the incident photons. Consider 
the the flat fold mirror to be at the origin, and the unit 
vector a points to the dewar, and b points to the posi- 
tion on the sky being observed. This is is illustrated in 
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FIG. 1. A pair of flat-fold mirrors with no net polarization. 
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FIG. 2. Fiducial geometry considered for calculating the relevant angles. Top left: the telescope beam from the dewar (direction 
a) to the sky (b), directed at the azimuthal polar coordinates (6 E ,(j>). Bottom left: the angle a through which the beam is 
reflected by the mirror. Right: E x and E y , the two polarization vectors on the sky. 



figure^ ol is the angle between a and b. In polar coordi- 
nates, a = (1, 90°, 180°) and b = (1, 6 e , 0) where 90 - 6 e 
and 4> are the the elevation and azimuth of the sky lo- 
cation. Converting into Cartesian coordinates, we have 
a = ( — 1,0,0) and b = (sin# e cos</>, sin# e sin cj), cos 8 e ). 
To determine a, we take the dot-product of a and b and 
get 



a.b 



cos a 



sin 8 e cos < 



(5) 



As the angles of incidence and reflection, 6^ and r are 
equal, we have 



l/2a = 1/2 arccos(— sin0 e cos (/>)) 



(0) 



We define E x to be parallel to the horizon, and E y to be 
"upwards" towards the zenith. We require the angle /3 by 
which E x and E y need to be rotated (around the line of 
sight) so that the rotated E x is parallel to the mirror. We 
define c, a unit vector that points in the same direction 
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as E x , so, in polar, coordinates, c = (l,90°,</> + 90°). 
aAb is parallel to the mirror. By taking the dot product 
of a A b with c, we can determine the angle j3 as 

(aAb).c = |(aAb)||c|cosj8. (7) 

After some algebra, 

r = arccos(cos 2 4>). (8) 

Thus we have the angle of incidence, 9i , required for the 
calculation above, and we have the angle, /3, so that the 
Rh and R± factors can be applied to the rotated E x and 



IV. RESULTS AND DISCUSSION 
A. Simulation Results 

We use CAMB to simulate a sky, observe it with our 
scanning strategy code, and follow the pipeline through 
to calculate power spectra. Firstly we test our pipeline 
by ignoring the effect of the fold mirror; resultant APS 
are shown in figure [3l In the region of interest {I ~100- 
1000) the output spectra match the input spectra well for 
TT and EE. The spectra from the masked pipeline maps 
for TT and EE are slightly supressed for larger £; 5% at 
at t ~ 900 when compared to the spectra made from the 
full sky maps. However, as the spectra from the masked 
full sky maps are also supressed, we conclude that this is 
due to the effect of the partial sky mask and fairly naive 
apodization procedure. Similarly, the BB spectra from 
the masked map is higher than the all-sky map, but as 
the pipeline BB spectra agrees with the masked all-sky 
map BB spectra, we conclude that it is the mask rather 
than the map-making procedure which is increasing the 
BB power. Although this level of analysis is not sufficient 
for a full B-mode experiment, this level of analysis is suf- 
ficient for investigating the effect of the flat-fold mirror. 

We then apply the effect of a flat fold mirror in the 
pipeline, using the resistivity of aluminium of 28.2 x 
10~ 9 ilm, and observations at 90GHz. A section of the 
time stream is shown in figure [4j both the input and out- 
put E% and the azimuth of the observed sky location. The 
correlation between induced polarization with observed 
sky azimuth is clear. 

The induced polarization is big, and so it is reasonable 
to remove it without requiring the entire data set. By 
applying the reverse of the observing algorthim, we can 
determine a "pre-mirror" E x and E y for a given mirror 
resitivity. We use £(-E^ — Ey), summed over the day of 
observation for a particular pixel, as a measure of the 
amount of polarization, and minimize it by varying the 
resistivity, using the Nelder-Mead simplex direct search 
13], implemented in Matlab. This time-line can then 
be re-inserted into the pipeline, maps made and spectra 
computed. The resultant APS are shown in figure [31 We 
see that in the range I ~ 100-1000, the TT, EE and BB 



spectra all match the pipeline output as before. In order 
to show problems more clearly, we also plot the difference 
between the spectra produced from the maps with and 
without the flat fold mirror. These are shown against 
I in figure El We find that the TT and EE spectra are 
reconstructed well, with less than 1% difference between I 
of 250 and 900. The BB spectra shows both fluctuation 
of ~ 1% and an overall offset of 2.5% over the same I 
range. 

B. Noise 

The shortest required observing time, for a general- 
ized CMB experiment, would be achieved using a photon 
noise limited telescope at a good site. For a instrument 
such as Clover, situated in the Atacama desert, the sys- 
tem NET is llhpK Vs. With a 100Hz data stream, this 
noise level overwhelms the induced polarization signal 
from T — > Q,U mixing. However, simple binning by a fac- 
tor of ten is sufficient to see the induced polarization sig- 
nal in the presence of white noise. If we add white noise 
at the level of 175 pK </s, with an p = 28.2 x 10" 9 fim, 
the best fitting resistivity returned with this procedure 
is 28.8 x 10" 9 fhn. 

Detector 1/f noise is assumed to be unimportant, as- 
suming a detector-based modulation strategy is imple- 
mented (e.g. using a rotating or stepped half-wave plate 
such as EBEX, or a fast detector switching such as 
QUIET). Atmospheric 1/f total power variations con- 
tributes strongly however as the effect of the reflection 
is a total power to polarization conversion. Thus the at- 
mospheric 1/f becomes a polarization 1/f term which 
does not get removed by any modulation strategy. A 
straightforward differencing strategy, for example differ- 
encing successive samples during azimuthal scanning, or 
differencing two adjacent detectors, to remove this vary- 
ing atmospheric contribution is likely to be the most ef- 
ficient method, but would reduce the effective observing 
time by a factor of two. 



C. Temperature variations 

The resistivity of aluminium varies with temperature, 
and so the induced polarization will vary with temper- 
ature. Given a resistivity at room temperature of po = 
28.2 x 10~ 9 f2m, the difference between the two polariza- 
tion reflection coefficients is 0.075%. We introduce a tem- 
perature coefficient of resistance, pr, of 3.8 x 10 _3 K _1 , 
with p = po(l + ptAT). Varying the temperature by 
20°C gives a difference of 0.078% between the two po- 
larization reflection coefficients, meaning that the T — > 
Q,U mixing varies by ~ 7%. For a reasonably large flat 
fold mirror - in this case the flat is ~ 2m diameter - the 
temperature variation is likely to be on the timescale of 
hours, driven largely by the ambient temperature vari- 
ations. Such variations should be removable by fitting 



■5 




FIG. 3. TT, EE, and BB spectra from various maps, dash-dot-dot-dot line is the input spectra, as simulated by CAMB. 
Dotted line is the all sky map, masked with the same mask used for the observed patches, solid is after our pipeline, but with 
no mirror, dashed is the spectra after having been observed with a flat fold mirror, dot-dash is the spectra after observed, but 
with the fitting algorthim (described in the text) applied. 



the resistivity to the data over different timescales on a 
pixel-by-pixel basis. We note that, apart from the pixel 
on the optical axis, each pixel will see a section of the 
mirror whose location varies in a predictable way with 
azimuth (and elevation) and so fitting with respect to 
mirror position might also be required. In this case, the 
error on the fit can be reduced by fitting multiple pixel 
data simulatiously to produce a mirror resistivity map. 



D. Oxide growth 

The polarizing effect of the fold mirror comes from 
its resistivity. As such, effects that alter the resistivity 
alter the polarizing effect. Aluminium grows an oxide 
layer very quickly, which then hydrates. The hydrated 
aluminium oxide can then be degraded chemically, or by 
physical damage to the surface. Such a surface is thin (5- 
15 nm) and so does not change the resistivity by a large 
amount. It is not unreasonable to assume that there will 
be changes to the surface during an extended observing 
campaign, both with oxide thickness and additional sur- 
face contamination. As with the temperature variation, 
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FIG. 4. Upper plot: time stream of the input and observed E x 
for a noiseless simulation, with a 28.2 x lCP 9 f2m fold mirror. 
Lower plot: time stream of the telescope azimuth. 
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FIG. 5. Fractional errors for the TT (dotted), EE (dashed) 
and BB (dot-dash) spectra. Calculated as the difference be- 
tween the no-mirror case and the fitted mirror case, divided 
by the no-mirror case. 

contamination would almost certainly be localised, and 
so fitting resistivity with respect to mirror position would 
probably be required. However such changes would be 
expected to be slow (timescales of days or more), or at 
least sudden and stable, and so more data is available to 



fit out the effect. As a result, we expect the resistivity 
variation from surface oxidation and contamination to be 
better characterised than the temperature variation. 

E. Polarized foregrounds 

The work so far has assumed that the light incident on 
the flat-fold mirror is polarized at a very low level, and 
that the majority of the polarized signal comes from the 
flat-fold mirror. This is the case when considering the 
CMB only, where the polarization signal has zero mean. 
These assumptions do not hold, however, when observ- 
ing an arbitrary patch of sky, where both synchrotron 
and dust from the galaxy is expected to be polarized at 
the ~ 1 - 10% level, with a total intensity of 10-100 
close to the Galactic Plane, but down to < l^iK at high 
Galactic latitudes. As such, we would recommend this 
method only be used in regions of known low foreground 
emission, where the resulting error in polarization recon- 
struction would be small. Over sufficiently large areas of 
sky the magnitude of the effect could be modelled, and 
in regions with relatively little large-scale variation in the 
foreground emission, the effect would be a constant off- 
set, similar to an additional absolute calibration error. 
The areas of sky observed by most B-mode experiments 
are selected to meet these requirements. 

In order to test this, we re-ran the simulation with 
some fiducial dust and synchrotron models added. After 
applying the same method as discussed above, the map- 
plane residuals are large; we see residuals of +20 to -40 
in the Q map, in comparison with the signal level of 
±20/xK across the map. We leave fully developing this 
method to work in the presence of significant polarized 
foregrounds for future work. 

V. CONCLUSIONS 

We have analyzed the effect of a flat-fold mirror on 
polarization measurements of the CMB. In the case of 
no foregrounds and no noise with no corrections, we find 
that the E and B spectra are heavily contaminated by the 
mixing of total power to polarization signal. However, as 
the CMB is expected to have a net polarization of zero, 
the effect can be fitted out. In the case of having a net 
polarization, such as in a real observation of the sky, the 
reconstructed maps contain significant residuals. 
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